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Superconductivity and spin-orbit (SO) interaction have been two separate emerging fields until very 
recently that the correlation between them seemed to be observed. However, previous experiments 
concerning SO coupling are performed far beyond the superconducting state and thus a direct 
demonstration of how SO coupling affects superconductivity remains elusive. Here we investigate the SO 
coupling in the critical region of superconducting transition on Al nanofilms, in which the strength of 
disorder and spin relaxation by SO coupling are changed by varying the film thickness. At temperatures T 
sufficiently above the superconducting critical temperature clear signature of SO coupling reveals itself 
in showing a magneto-resistivity peak. When T < T^, the resistivity peak can still be observed; however, its 
line-shape is now affected by the onset of the quasi two-dimensional superconductivity. By studying such 
magneto -resistivity peaks under different strength of spin relaxation, we highlight the important effects of 
SO interaction on superconductivity. 

Spin-orbit (SO) interaction, which couples the electron orbital motion to its spin, has been at the center of 
intensive research efforts in the field of spintronics^"^. This coupling between orbital motion and spin 
orientation leads to several interesting physical phenomena such as the spin Hall effects^ ^ spin- orbit gap^'^ 
and spin relaxation^'^, and may open an avenue to achieve all -electrical control of spin degree of freedom^^. In 
addition, SO interaction provides an essential ingredient in creating topological insulators and Majorana 
fermions^^. Furthermore, recently SO coupling was found to have a significant impact on quasi two-dimensional 
(2D) superconductivity. One of the most fascinating examples is the promotion of superconductivity with 
enhancing the strength of SO interactional In the 2D electron gas at the interface between LaAlOs and SrTiOa 
exhibiting superconducting properties at low temperatures T, large Rashba SO coupling arising from interfacial 
breaking of inversion symmetry is shown to have pronounced effects in stabilizing the superconducting state^^'^^. 
Moreover, due to strong SO coupling at the interface, ferromagnetism, which is usually thought to be incom- 
patible with conventional superconductivity, can now coexist with superconductivity Interestingly, SO coup- 
ling also appears relevant to the observed increase of superconducting critical temperature by an in-plane 
magnetic field in Pb thin films^^. Within the framework of the Bardeen-Cooper-Schrieffer (BCS) theory^^ which 
attributes superconductivity to pairs of electrons with opposite spin projections (Cooper pairs), any mechanism 
that causes spin alignment (like those induced by paramagnetic impurities and magnetic fields) is not predicted to 
enhance superconductivity. Response of superconductivity to magnetic perturbations can be fundamentally 
different in the presence of SO coupling. 

In a quasi 2D superconductor with SO coupling, the application of a magnetic field perpendicular to the plane 
of carrier motion may not only lead to suppression of weak antilocalization (WAL) induced by SO coupling^^"^^, 
but also to orbit de-pairing of Cooper pairs, destroying the zero -resistivity superconducting state. Magneto- 
transport studies could thereby allow one to obtain fundamental information regarding the interplay between SO 
coupling and superconductivity. Here we present such measurements on molecular-beam epitaxy (MBE) -grown 
Al nanofilms with thicknesses of 3 nm, 6 nm, and 12 nm in which the rate of spin relaxation by SO coupling 
depends on the nanofilm thickness. To date, most measurements concerning SO coupling are performed far 
beyond the superconducting state, making a direct demonstration of how SO coupling affects superconductivity 
unavailable. In our experiments, when T > Tc, a WAL peak in the magneto -resistivity due to the presence of SO 
coupling is observed in all the nanofilms. By probing the temperature evolution of this resistivity peak toward the 
superconducting state, we have demonstrated that WAL effects are progressively suppressed as a result of the 
onset of superconductivity. In the superconducting state, the peak shape is clearly affected by Cooper paring. 
Moreover, an unusual positive contribution to the magneto -resistivity (unexpected within the BCS theory and 
WAL effects) appears near the resistivity peak for the 3-nm-thick Al film in which the relaxation of spin is fastest 
among the three studied films. Our experimental work therefore highlights the important effects of SO coupling 
on superconductivity. 
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Figure 1 | (a) Optical image of the 6-nm-thick Al film, (b) Longitudinal resistivity Pxx as a function of T. (c) Longitudinal and Hall resistivities Pxx and 
as a function of 5 up to 11.7 T at T = 4.4 K. p^{B) at various T. From top to bottom at 5 = 2 T: (d) T = 24 K, 22 K, 20 K, 18 K, and 16 K; (e) T = 8 K, 
10 K, 12 K, and 14 K; (f) T = 3 K, 4 K, and 6 K. The red curves correspond to the best fits to Eq. (1). The inset of (f) presents the fit considering 
Maki-Thompson superconducting fluctuations to the data at T = 3 K. 



Results 

SO coupling in the normal state. Figure 1(a) shows the optical 
image of the studied 6-nm-thick device and Fig. 1(b) shows the 
obtained longitudinal resistivity Pxx over a wide range of 
temperature (5 K < T < 30 K). For T > 18 K, Pxx increases with 
increasing T, characteristics of a metal film subject to electron - 
phonon scattering. For T < 18 K, quantum corrections primarily 
due to WAL and superconducting fluctuations would then instead 
govern the T dependence of resistivity as the effects of electron - 
phonon scattering is suppressed. Figure 1(c) shows the magneto- 
resistivity measurements for 0 < 5 < 11.7 T at T = 4.4 K. We can 
observe that Pxx is much larger than the Hall resistivity Pxy over the 
whole field range. The smallness of Pxy is due to a large carrier density 
in metal systems. Therefore one is able to obtain the longitudinal 
conductivity of ^ II according to = p^/{p^^ + Pxy')- 

Figures l(d)-(f) show Pxx(^) at various T. As presented in Fig. 1(e), 
for 8 K < T < 14 K, signatures of WAL effects, that is, positive 
magneto -resistivity (MR) around B = 0 followed by negative MR 
with increasing B, are observed, clearly indicating the importance of 
SO coupling in the measured Al nanofilm with thickness of d = 
6 nm. For T > 16 K, weak localization (WL) instead of WAL dom- 
inates the transport process, inferred from the monotonically decreas- 
ing trend of Pxx with the applied B shown in Fig. 1(d). Suppression of 
WAL at elevated T is expected to occur when no significant spin- 
dependent phases can be accumulated. In metals, the dominant 
mechanism of SO coupling is believed to arise from the potential of 
lattice ions such that the spin orientation is randomized following the 
Elliot- Yafet (EY) mechanism which states that spin flips during every 
scattering event^'^^'^^. According to the WAL theory developed by 
Hikami, Larkin, and Nagaoka (HLN) (ref. 26), which assumes that 
the EY mechanism is responsible for spin relaxation, the MR in the 
2D diffusive regime under 1/pxx has the form 

^2 



where crj^ is the zero -field conductivity independent of B, ^(x) is the 

4 

digamma function, Bi = Bq + Bso^ B2 = B^^ + - (Bso)y and B^ = B^. 

Here the characteristic fields Bq, Bso and 5^ are related to the elastic 
scattering length Iq, spin- orbit relaxation length /so» and dephasing 
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curves in Figs. l(d)-(f) denote the best fits to the data using Eq. (1). 
The excellent agreement between experiments and HLN theory 
shown in Figs. 1(d) and 1(e) for 8 K < T < 24 K suggests that the 
EY process is indeed the dominant mechanism for spin relaxation. 
However, it is found in Fig. 1(f) that HLN theory fails to describe the 
MR behavior for T < 6 K, which is ascribed to the influence of 
fluctuating Cooper pairing (superconducting fluctuations). The inset 
of Fig. 1(f) presents the fitting result considering SO interaction and 
Maki-Thompson superconducting fiuctuations. A good agreement 
can likewise be observed. However such a method is not applicable 
as the superconducting state is well developed. 

SO coupling in the superconducting state. In order to study the 
interplay between SO coupling and superconducting fluctuations 
around 5 = 0, we limit the fitting range to high magnetic fields, in 
which superconductivity is completely suppressed, to obtain the SO 
relaxation and dephasing length. Here we choose the field B^ at 
which Pxx(^) reaches its maximum as the lower bound to fit our 
data to Eq. (1). The underlying reason for this choice is that 
superconductivity cannot give rise to negative MR. On the other 
hand, the upper limit of this fit is chosen at 5 = 2 T to avoid the 
influence of positive MR which occurs at high fields due to orbital 
motion as shown in Fig. 1 (c). The new fits to the data are shown as the 
red curves in Fig. 2(a) for 3 K < T < 6 K and further in Fig. 2(b) for 
0.328 K < T < 1.980 K, where remarkable agreement between 
theory and experiments is found. Three important length scales Iq, 
Iso, and can thereby be determined and the results are plotted 
against T in Figs. 3(a) and 3(b). We observe that /so exceeds for 
T > 14 K, indicating that phase coherence has lost before spin- 
dependent phase is accumulated for T > 14 K. Consistently, WAL 
effects are barely observable for T > 14 K in Fig. 1(d). As shown in 
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Figure 2 | pJ^B) at selected temperatures Tfor the 6-nm-thick film. From top to bottom at 5 = 2 T: (a) T = 3 K, 4 K, and 6 K; (b) T = 0.328 K, 0.914 K, 
1.568 K, and 1.980 K. The red curves represent the best fits to Eq. (1) over the limited field range (5p < 5 < 2 T) and the blue curves denote the 
extrapolations of the fits to low B regime using the extracted parameters. The black and blue arrows indicate the experimentally obtained and theoretically 
predicted peak positions. To emphasize the zero-resistivity superconducting state around 5 = 0 when T < 1.980 K, the fijU-scale view at the 
corresponding Tare shown in the inset of (b). (c) Pxx(T) at various B around the superconducting transition. 



the inset of Fig. 2(b), the zero -resistivity superconducting state 
appears when T is low enough. The characterization of the 6-nm- 
thick superconducting nanofilm which yields important physical 
quantities can be found in the supplementary information. To 
further study the effects of Cooper pairing, we plot the theoretical 
curves given by Eq. (1) around B = 0 with the parameters determined 
from the fitting procedure at high B as described above. These 
simulation curves, which are the extrapolation of the fitting curves. 
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Figure 3 | Dephasing length l^^ spin-orbit relaxation length /so> and elastic 
scattering length Iq as a function of T ranging from (a) 8 K to 24 K and 
(b) from 0.328 to 6 K. (c) The difference between the measured and 
theoretically predicted resistivity peak position (ppandpp ^^at B = B^ 
and5p respectively) as a function of T. 



are represented by the blue ones in Fig. 2(a). The expected infiuence 
of superconductivity on WAL can then be tracked as a function of T 
even down to 0.328 K, well inside the superconducting regime. By 
plotting Pxx(^ at different B in Fig. 2(c), we observe that there is a 
crossover between superconducting and insulating states with B. 

By comparing the experimental data (where both SO and super- 
conductivity exist) with the simulated curves (where only SO is con- 
sidered) as shown in Figs 2(a) and 2(b), we observe that the position 
of resistivity peak, characterized by strong SO coupling in the normal 
state, is shifted to a higher magnetic field compared with the expected 
position when T < 6 K. Figure 3(c) then collects the difference 
between the measured and theoretical peak positions in their res- 
istivity and magnetic field values at various T covering both the 
normal and superconducting states. This difference becomes signifi- 
cant at low T, suggestive of suppression of WAL effects by spin- 
singlet Cooper pairing. In the normal state, positive MR around B 
= Ois due to SO coupling. In contrast, in the superconducting state, 
orbit de-pairing by magnetic field, leading to notable positive MR, 
plays a decisive role in the low-field magneto-transport properties 
since SO coupling has limited effects on carriers with total spin zero 
such as bosons. However, in the transition regime from normal to 
superconducting state, the orbital motion of spin-half electrons, 
which are fermions, still couples to their spin orientations due to 
SO interaction. 

As shown in Fig. 2 (b), in the normal state indicated by the red 
fitted curves, Pxx decreases with increasing temperature, character- 
istic of a weak insulator. When B is below 5p, positive magneto- 
resistivity is observed and insulating behavior persists before the 
resistivity of the device reduces dramatically due to superconduct- 
ivity. Since the Pxx maxima predicted by the standard HEN model 
(indicated by the blue curves) occur at a substantially lower B than 
those observed in experiments, Pxx for ^ < cannot be solely 
described by spin-orbit coupling in the normal state and thus must 
be related to superconductivity. Such results are reminiscent of the 
intriguing infiuence of disorder on super conductivity^^"^^. In a dis- 
ordered superconductor subject to a strong perpendicular magnetic 
field, insulating behavior can occur. For a moderate amount of dis- 
order in the system, the observed insulating behavior at high fields is 
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Figure 4 | (a) & (b) Magneto-resistivity measurements Pxx(5) at different temperatures Tfor the 3-nm-thick film. The red curves correspond to 
the best fits to Eq. ( 1 ). (c) - (f) Zoom-in of (b) at each separate T. The insets show Pxx (5) around Bp indicated by the arrows. The red curves are the best fits 
to Eq. (1) over the Umited field range (Bp < B < 2 T) and the blue curves denote the extrapolations of the fits to low B regime using the extracted 
parameters. 



generally ascribed to localization of unpaired electrons. However, in 
a strongly disordered superconductor, increasing magnetic field does 
not always quench superconductivity^^"^^. Instead, thanks to the 
formation of superconducting islands (SCIs), with the presence of 
strong disorder, an applied out-of-plane magnetic field, which sup- 
presses the correlations between different superconducting islands, is 
not detrimental to superconductivity drastically. The route toward 
insulating behavior with increasing magnetic field thereby proceeds 
through development of well -separated superconducting islands fol- 
lowed by localization of unpaired electrons, which is unexpected 
within the conventional BCS theory. Such an unusual superconduc- 
tor-insulator transition (SIT) has already been observed in a wide 
variety of amorphous superconducting thin films including InxOy 
(ref. 27), TiN (ref. 29), and Bi (ref. 31). Strong disorder, which is 
believed to be detrimental to superconductivity, plays a crucial role in 
maintaining superconductivity in such systems. In the presence of 
strong disorder, normal -state transport should be governed by 
activation or hopping process, giving rise to an exponential T 
dependence of resistivity^^'^^'^^'^^. In our MBE-grown nanofilms, dis- 
order may not be strong enough to support the formation of SCIs, 
which can be understood from the observed weak T dependence of 
Pxx at B = 2 T shown in Fig. 2(b). Moreover, the obtained Pxx over 

h 

the whole measurement range is far below — - ~ 13 kQ which is the 

critical resistivity value expected for the onset of strong localization. 

Further magneto -transport measurements were performed on the 
Al nanofilms with thickness of 3 nm and 12 nm, which are expected 
to have respectively faster and slower spin relaxation compared with 
that of 6 nm. Figures 4(a) and 4(b) show the measured Pxx(^) at 
various T for the 3-nm-thick film. We observe that positive MR 
around B = 0 due to WAL, can be well fitted to Eq. (1) for T > 
8 K. For T < 8 K, influence of superconductivity on the magneto - 
resistivity becomes increasingly important with decreasing T. For 
clarity, we plot Pxx(^) at T = 5 K, 4.2 K, 3 K, and 2 K separately 
in Figs. 4(c) -(f). The experimental data are in good agreement with 
WAL theory described by Eq. (1) when B > Bp, the field where the 
resistivity peak occurs. Following the same analysis procedure done 
for the 6-nm-thick film, we extrapolate the high-field fitting curves to 
lower fields by directly substituting the extracted fitting parameters 
into Eq. (1). Interestingly, as seen by the comparison between mea- 
sured and simulated curves for B < Bp, superconductivity does not 



always cause the measured resistivity to be lower than that predicted 
by Eq. (1) for a given B. For example, figure 5(c) shows that at T = 
5 K the magneto-resistivity near Bp is higher than its expected value 
without taking superconducting effects into account given by Eq. (1). 
With decreasing T, this unusual positive resistivity contribution 
becomes gradually smaller and eventually disappears, which is 
observed clearly in the insets of Figs. 4(c) - (f), the zoom-in of the 
region around Bp. At lower T and B, the resistivity value is always 
lower than the prediction of Eq. (1) due to stronger Cooper-pair 
correlations. Such results indicate that there will be a mechanism 
relevant to Cooper pairing which makes the resistivity higher than 
that predicted by Eq. (1). 

The enlargement of SCI size with decreasing B is thought to be one 
of the mechanism which can provide an extra resistivity contri- 
bution. With a decrease in B, the enlarged SCIs may further block 
the electron transport. As observed in Fig. 4(b), even for the 3-nm- 
thick film, the T dependence of Pxx around Bp is still very weak and 
h 

Pxx (B) ^ — r —13 kQ, suggesting that disorder does not play a key 

role in our observed resistivity anomaly. It is noted that the magneto - 
transport measurements followed by the same analysis procedures 
have also been done for T < 2 K. All of them share the same features 
with that at T = 2 K shown in Fig. 4(f), where Cooper pairing 
appears to cause the measured resistivity to be lower than that pre- 
dicted by Eq. (1) over the whole B range below Bp. At T < 2 K, since 
superconducting correlations are strong, this unusual resistivity con- 
tribution is completely suppressed. Finally figures 5(a) and 5(b) show 
the magneto -transport data measured on the 12-nm-thick film. The 
resistivity peak well described by Eq. (1) in the normal state is com- 
pletely suppressed at low temperatures. In the 12-nm-thick film, 
strong superconducting correlation completely suppresses the influ- 
ence of SO coupling. 

Discussion 

The physical quantities of superconductivity and SO coupling 
obtained from the fits for the Al nanofilms of various thicknesses 
are summarized in Table 1. We note that ^ is longer than the film 
thickness in all our samples. Therefore we always have quasi-2D 
superconductivity. We can clearly see in Table 1 that Tso (— ho^ 
D) (ref 38) decreases with decreasing the film thickness, indicating 
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Figure 5 | (a) & (b) Magneto -resistivity measurements Pxx(5) on the 12-nm-thick film at different temperatures T. The red curves correspond to the best 
fits to Eq. (1). 



that the strength of spin relaxation is stronger in a thinner film. 
Correspondingly, the amount of disorder increases with decreasing 
the film thickness, which can be understood from the comparison of 
the zero-field resistivity at T = 8 K well above and from the 
temperature dependence of resistivity p(r, 5 = 0) for T > 8 K (see 
the Supplementary Fig. SI). We can see clearly that the insulating 
behavior becomes progressively stronger with reducing the film 
thickness. It is known that two extensively studied effects of spin 
relaxation are the Elliot- Yafet and Dyakonov-Perel (DP) mechan- 
isms^'^^'^^'^°. In the EY mechanism, spin relaxation is relevant to how 
the spin changes its orientation during a scattering process. This is 
possible as the electron wave functions associated with a given spin 
acquire an admixture of the opposite-spin state due to SO coupling 
induced by lattice ions. In contrast, the DP mechanism is related to 
the Dresselhaus SO coupling caused by bulk inversion asymmetry 
and the Rashba SO coupling caused by structural inversion asym- 
metry. In this case, a charged particle moving in an electric field 
(originating from the lack of inversion symmetry) experiences an 
effective magnetic field, which couples to its spin. The spin relaxation 
is thereby relevant to spin precession around this effective magnetic 
field between scattering events. Therefore, increasing disorder may 
enhance the rate of spin relaxation govern by EY mechanism. 
Moreover, as shown in the Supplementary Fig. S2, the model 
developed by lordanskii, Lyanda-Geller and Pikus (ILP) (ref 41), 
in which the DP mechanism is included, cannot yield good fits to 
the measured MR data as compared to the HLN model especially in 
the 3-nm-thick film with the strongest disorder in our experiments. 
These results indicate that EY mechanism dominates over DP one 
with decreasing the nanofilm thickness. 

The scale of energy variation involved due to spin relaxation can be 
evaluated utilizing the uncertainty relation Sso — h/Tso (ref- 14). By 
calculating the ratio of Sso to 2z/, which is found to be 0.37, 0.77, and 
0.85 for the films of d = 12 nm, 6 nm, and 3 nm, respectively, the 



growing importance in the correlation between spin relaxation and 
superconductivity with decreasing the film thickness is thereby 
revealed. Correspondingly, /so approaches progressively to ^. This 
numerical analysis suggests that the unexpected resistivity contri- 
bution around Bp observed in the 3-nm-thick film is indeed related 
to the combination effects of SO coupling which causes spin relaxa- 
tion and superconductivity which pairs electron spins. It is well- 
established that the incoming electron with spin-up takes another 
electron with spin-down to enter the BCS superconductor as a 
Cooper pair with zero spin. The attractive interaction that pairs 
the electrons to form superconducting state can be mediated by 
lattice vibrations as described by the BCS theory. Since the SO coup- 
ling involved here is provided by the potential of charged atoms, in a 
material preserving ideal crystalline structure, spin orientation 
relaxes spatially homogeneously. However, in the presence of dis- 
order this process will then become inhomogeneous. Our experi- 
mental findings thus provide an systematic way to probe the 
influence of SO coupling on superconductivity by studying the 
evolution of magneto -resistivity peak in the superconducting state. 

In conclusion, magneto -transport in the Al nanofilms of different 
thicknesses (3 nm, 6 nm, and 12 nm) was investigated using weak 
antilocalization analysis. For all the films in the normal state, we 
observe clear signatures of WAL, that is, positive magneto-resistivity 
around B = 0 followed by negative one with increasing 5, indicating 
the importance of spin-orbit coupling. Moreover, all the measured 
MR in the normal state agrees well with the calculation work done by 
Hikami, Larkin, and Nagaoka. With decreasing temperature, which 
causes the onset of superconductivity, we show that the MR peak due 
to SO coupling in the normal state would now be govern by the 
interplay of SO coupling and superconductivity. For the 6-nm-thick 
film. Cooper pairing reduces the resistivity from the value predicted 
by HLN theory and thereby shifts the MR peak to a higher B. More 
interestingly, we find that there is an extra resistivity contribution 



Table 1 | Summary of the physical quantities of superconductivity and spin-orbit coupling 
Superconductivity Spin-orbit interaction 



Disorder strength Coupling strength 



2A 
meV 



Hc(0) 
T 



D 

10-^ mVs 



10- 



eso 
meV 



pJT=8KB- 



0) 



o/2A 



3 1.509 0.46 0.339 31.18 2.44 51.30 1.08 0.39 1099 0.85 1.65 

6 1.577 0.48 0.424 27.88 2.04 48.12 1.13 0.37 154 0.77 1.73 

12 1.358 0.41 0.139 48.69 5.36 120.90 2.72 0.15 22 0.37 2.48 



°/so shown here is obtained by averaging the extracted parameters below 4.5 K which are normally T-independent. 
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close to the MR peak occurring accompanied by the appearance of 
superconductivity in the 3 -nm- thick film having faster spin relaxa- 
tion. For the 12-nm-thick film, since superconducting correlation is 
rather strong and spin relaxation is comparably slower, MR peak is 
completely suppressed at low T. Different behavior in the films of 
different thicknesses results from different strength of spin relaxation 
and superconductivity. Spin-orbit coupling, which makes the spin 
degree of freedom respond to its orbital characteristics and causes 
spin relaxation, has thereby been demonstrated to have important 
effects on superconductivity. Our experimental work sheds new light 
on the understanding of the role of SO coupling in a superconductor. 

Methods 

All samples were grown in situ using a Varian Gen II solid- source molecular beam 
epitaxy (MBE) system. For an Al film whose thickness is 30-100 nm, an undoped 
200-nm-thick GaAs buffer layers is grown on semi-insulating (100) GaAs substrate as 
a template. We find that the thickness of the GaAs buffer layer needs to be increased to 
300 nm to achieve a flat enough template for depositing an aluminum nanofilm 
which is thinner than 15 nm. Before Al deposition, the buffer layer on top of the semi- 
insulating (100) GaAs substrate is kept in the ultra-high vacuum to prevent the GaAs 
surface from oxidation. After the growth chamber is pumped down to 3 X 10"^° Torr 
to eliminate residual arsenic vapor pressure, an Al nanofilm is grown at ~0°C at a 
slow rate of 0.18 |j,m/h (0.05 nm/s) (ref 42). We note that once the film thickness is 
below 2 nm, the Al nanofilm does not conduct at room temperature, suggesting that 
an ultrathin Al film gets oxidized and become AI2O3 which appears to be an insulator. 
For an Al nanofilm whose as-grown thickness > 3 nm, the top AI2O3 film serves as a 
protective layer so that the underneath Al nanofilm does not get oxidized during 
processing. All samples were processed into 50-|im-wide Hall bars using conven- 
tional photo -lithography and lift-off technique. The mesa is defined by 2% tetra- 
methylammonium hydroxide (TMAH) etching for 20 seconds for pattern 
transferring from photoresist to aluminum film. It is found that such a slow etch rate 
is crucial for high-yield, manageable and reliable sample preparation. A 30 nm/ 
300 nm Ti/Au layer is deposited using electron-gun evaporation as the contact 
electrodes. Four-terminal magneto-transport measurements were carried out using 
dc constant-current sources in a top-loading He^ cryostat equipped with a super- 
conducting magnet. Since the device was patterned into Hall bar geometry one can 
obtain both longitudinal and Hall resistances and R^. In Fig. 3(a), we know that 
the mean free path Iq is larger than the thickness (~6 nm) of the film being measured, 
indicating that the nanofilm in the present study can be treated as a quasi-two- 
dimensional system. Hence we can define longitudinal and Hall resistivities and 
Pxy by Pxx = Rxx^/l and = R^, where w is the width of the Hall bar and / is the 
distance between the two voltage probes used for the measurement of R^x- The same 
results were obtained on 3-nm and 12-nm-thick films. Source-drain current /sd on 
the order of 10 |J,A was injected to improve the signal-to-noise ratio. For I^u < 1 M^A, 
the weak peak structures were hardly observed. Although large current may cause 
heating effects, which increases the electron temperature, our main findings 
regarding the coupling of spin-orbit coupling and superconductivity would not be 
influenced since /sd is not large enough to completely destroy the superconducting 
state. 
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